5 Symmetry 

5.1 Rotation Operators 

Consider a scalar field f(r). 

Example 1: The temperature at point r 

After a rotation of the coordinate system, the same field is described by 
the new function f'(r'), with 



Here the operator R acts on the coordinates. Introduce now the operator 
Pr, which acts on the function / in Hilbert space, making it follow the 
rotation of the coordinate system: 



Active interpretation: In the left-hand equation, f'(r') = Pnf(r'), the 
operator Pr changes the function / into /'. 

Passive interpretation: In the right-hand equation, rewritten as f(r') = 
Pr 1 fin), the operator does not change the function but acts on the coordi- 
nates r, changing them to r' 

Example 2: Let f{<p) = e im< ^ and R(f> = <p' = (f) — a, a positive rotation 
through a: 



r' = Rr 



(5.1) 



/'(/) = P R f(t) = f(r) 



(5.2) 





\ 



\ 



\ 



\ 



\ 




Figure 5.1: Positive rotation through a 
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Then 

e im(0'+a) _ e ima ^nup' _ Q \m<t> ^g g-j 

/'(</>') ^fl /(0O /(<« 

5.1.1 The Euler Angles 

Consider a right-handed coordinate system Oxyz and positive rotations (like 
in fig. (5.1), R z pushing a right-handed screw towards +z). Introduce the 
rotations 

(i) Rotation R z (a) through a around O z (O y — > O y /) 

(ii) Rotation R y >(f3) through (3 around O y i (O z — > O z ») 

(hi) Rotation R z ii("f) through 7 around O z »(O y i = O y n — > Oy'») 

Then eq. (5.3) (replacing m by the operator j 2 ) 

/V) = e^'/d:') (5.4) 
/"(/') = e^V7'(rO 

= e i(3j y"e iaj *"f(r") (5.5) 
/'"(/") = e i7J '«"7"(r") 

= e i7i z »' e ift H »' e ia ^'" /"(/") (5.6) 

in the active interpretation using temporary coordinates O z , O y . 
The Pr operator of the type (5.2)-(5.3) becomes 

P R (a/37) = e i7iz e i/3 ^ e iQ ^ (5.7) 

With this operator, the same coordinates must be used on both sides. 

If fixed rotation axes, Oxyz, are used, one should replace (j^a) by (a/37) 
in the operator Pr. 



Rotations 


Axes 


Operations 






i.{z) 


2.(j/) 3.(z) 


Active 


Fixed 


7 


a 


Active 


Temporary 


a 


7 


Passive 


Fixed 


—a 


-/3 -7 


Passive 


Temporary 


-7 


-/? -a 



5.1.2 Rotation of Spherical Harmonics. 

In the active, temporary picture, the rotated function becomes 

Z'{f<f/) = e- i ^'e ift V e iQ ^'Yi rn (i? / 0') (5.8) 
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In the passive, temporary representation 

Y, ro (0V) = e- [a ^e-^ye-^Y lm m 

= ^D^, m (a/3 7 )Y^W) (5-9) 

m' 

Here the rotation matrices ( Wigner matrices) 

D l m , m = (Im'le-^'e-Wye-W'llm) (5.10) 
The temporary rotations l.a, 2./3, 3.7 bring Oxyz to Ox'y'z'. 

5.1.3 Rotation of |jm) Functions 

Because the functions, with half-integer j, obey the same commutation 
rules as those with integer j, we may simply substitute I for j: 

P R\jm) =E D !»'m(«WH ( 5 - U ) 
m' 

5.1.4 The 2-to-l Homomorphism from SU{2) to 50(3) 

The special unitary group, SU(2), consists of all unitary (U^ = U _1 ), 
unimodular (det U = 1) 2 x 2 matrices. We may regard it as the group of 
all three-dimensional rotations, spanned by the spin vectors 



(5.12) 



m = a = \--) 
m = P = \\-\) 

A Pauli matrix Oi transforms under the rotation U to 

x,y,z 

Uct;U 1 ^o-jAj, (5.13) 

3 

The matrix A is real, orthogonal (A^ A = 1) and unimodular and belongs 
hence to the special orthogonal group SO (3) (consisting of 3 x 3 matrices), 
isomorphic to the group R(3) of all three-dimensional rotations. 

RmjR* = ^ rrijAjj (5-14) 

3 

where the vectors rrij form a cartesian basis. 

Consider now two elements, U and V, of SU{2) and two elements, A 
and B, of SO{3), -► R(A) and 12(B). A comparison of (5.13) and (5.14) 
yields the mapping 

U^A (5.15) 
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Letting V -> B, (5.13) => 

(UVj^^Ut) = Uo-jBj-j (5.16) 

Hence (5.15) is a homomorphism. In particular, as 

U^Ut = (_U)(7 i (-U) t = o-j-Aji =► ±U -»• A, (5.17) 

it is a 2-to-l homomorphism, exactly two elements of SU(2) is mapped to 
exactly one element of 50(3). 

As the angular momentum j = a/2 for SU(2), the matrix operator for 
a rotation through <j> around n becomes (recall (5.3)): 

U(n, <f>) = S^ 12 = cos ^ + i{a ■ n) sin ^ (5.18) 

Hence 

U(n,2vr) = -1 (5.19) 
U(n,4vr) = 1 (5.20) 

5.2 Double Groups 

5.2.1 Non-relativistic Case with Spin 

Let G be the non-relativistic or simple group. Its elements, like the pure 
rotations, B4, commute with H, 

[H,Ri}=0 (5.21) 

The spin space is entirely decoupled from coordinate space. The total sym- 
metry group is the direct product 

Gfuii = G x SU{2) (5.22) 

For a Coulomb potential, Gf u n is larger still (Fock (1936)). 

5.2.2 Relativistic Case 

Now the space and spin variables are coupled. Only those Uj G SU(2) which 
correspond to an R« (by (5.15)) can be included in 

G* = {±RiUi} c G x SU{2) (5.23) 

The two signs correspond to those of (5.17). Let K be the set of all such 

K = {±LbjU -» RJ (5.24) 
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K and G* are isomorphic, K = G* . The relations between the four groups 
are seen in figure (5.2). 




Figure 5.2: Group relations 



5.2.3 Improper Rotations 



The improper groups, G', 



are either direct products 



a = g x d 



(5.25) 



where Cj = {E, 1} contains the inversion, /, or isomorphic to G. In the first 
case G' contains /. 

Example 3: 

T h = TxC t = {E, 3C 2 , 4C 3 , 3C^}V{I, 3a, 4S 6 \ AS 6 } (5.26) 
5.2.4 The Group 0(3) 

All improper groups, G', are subgroups of the group 0(3), which is obtained 
from SO (3) by adding the improper rotation with det A = ±1. 




If the inversion is included, it commutes with all other elements, 



V 



1 



being a multiple of I 



1 



V V 
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5.2.5 The Double Group is a Symmetry Group of the Dirac Equa- 
tion. 

Consider 

/i D = ca ■ p + (3mc 2 + V (5.27) 

with 

a,J°A. pJ 1 °) (5.28) 
and define the transformation matrix for the i th rotation as 

UP = f Ui ° ) (5.29) 

Consider the group 

G* D = {±RiUf } ^ G* (5.30) 
This rotation transforms the kinetic energy to 

(RTjP)a.p(R,UP)t= f ° UCTfeUt )R Pfe Rt (5.31) 

\Ua fe Ut J 

= aip m Ai k A mk = a m p m = a-p 

Thus the kinetic energy remains invariant under the double group opera- 
tions. 

5.2.6 The Element E 

Bethe (1929) arrived at the double-group concept by including in the group 
the new element E, corresponding to rotations through 2tt, for which 

E 2 = E (5.32) 

Example 4: 

C 2 = {E, C 2 }, C* 2 = {E, C 2 , E, C 2 = EC 2 } (5.33) 

In general 

A = EA (5.34) 
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5.2.7 Elements of Double Groups 

In the multiplication table of a double group the new, barred, operations 
are mixed with the others. For rotations around the same axis 




Figure 5.3: Multiplication of rotations in a double group. From Koster el 
al. (1963). 

The properties of the inversion are 

// = E, 11 = E, II = E, IK, = Sj, Itij = Sj (5.36) 
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As 

we get 



Also, 



a h = IC 2 (5.37) 

a h C 2 =lCl = IE = I + l\ (5.38) 

aa = {IC 2 ){IC 2 ) = C\ = E (5.39) 

aa = aaE = EE = E =>- (5.40) 

a- 1 =a\ (5.41) 

{E,a}® {E,E} = {E,a,E,a^} (5.42) 

aa' 1 = E (5.43) 

aE = a (5.44) 

Ea' 1 = Ea = EEa = a (5.45) 



For a rotation C n , 



(% = E, C 2 n n = E^ (5.46) 

EC n = C'nC n = C n C'n = C n E (5.47) 

Furthermore, 

f R ? , for < Rj < 7T 

C 2 R» = < 3 5.35 

1 Rj , for — 7r < Rj < 1 J 



[ Sj, for < Rj < 7r 
a h R = IC 2 Ri = < 5.48 
\Si, for-7r<R<0 V ' 

5.2.8 Irreducible Representations ("irreps") 

The order of the double group, #G* = 2(#G) = 251. Hence its irreps must 
satisfy 

r 

J>f = 2 5 (5.49) 
i=l 

r being the number of irreps for G* . 
Example 5: 
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Both C 2 h and C 2v have g = 4. They satisfy (5.49) by 

C 2h ■ (I 2 + l 2 + 1 2 +1 2 )+(1 2 + l 2 + l 2 + 1 2 )= 8 (5.50) 
C 2v : (l 2 + l 2 + 1 2 +1 2 )+ 2 2 =8 (5.51) 

ri ... r 4 r 5 

Thus all relativistic molecular orbitals of C 2v belong to the same irrep, T5. 
Note that, for an even number of electrons, the total electronic symmetry 

still belongs to T\ r4, J being an integer (for higher symmetries, when 

applicable) . 

The reps of G are also reps of G* 



We have a 2-to-l mapping from G* to G and a n-to-1 mapping from G 
to reps {T}. Thus we have a 2n-to-l mapping from G* to a single matrix of 
the representation T. 

The original irreps of G have basis functions, \jm), with integer j. For 
their characters 

X (i) (A) = x (i) (A) (5.52) 

The additional irreps of G* have basis functions \jm) with half-integer j. 
For them 

X W (A) = - X «(A) (5.53) 
Here i is the irrep in question. 

5.2.9 Classes 

The number of classes equal the number of irreps, 

k = r (5.54) 

The character of the rep, spanned by the functions \jm), m = —j, . . . , j, 
becomes 

<j), x sin(j + \)uj 

X u> {uj) = — (5.55) 

sin 

For uj = 7r and half-integer j, x^'H 71- ) = 0- Therefore x^K 7r + 2vr) = 
— X^H 71 ") = 0) as wen - Thus x(C*2) = X(C*2) C2 and C*2 may or may 
not belong to the same class. 

For n > 2, one always has x(C n ) = —x(C n ) and here C n and C n belong 
to different classes. 
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5.2.10 Theorem of Opechowski 

Opechowski (1940): C n and C n belong to to different classes if n > 2. For 
n = 2, C2 and C2 belong to the same class if and only if the group contains 
another C2, perpendicular to the C2 considered, or if it contains a symmetry 
plane, a, containing the C2 considered. 

Similarly, S3 and S2 belong to different classes except if there is a a, 
containing the C2 axis, or another C2 perpendicular to it. 

The reflections a and a belong to different classes, except if there exists 
another a, perpendicular to it, or a C2 in the plane a. 



Example 6: 

Character table and basis functions for the group Dg h : 



D 3 h 


E 


E 




2C 3 


2C 3 


2S 3 


2S 3 


3C£ 
3C^ 


3<T V 
3(T V 


Time 
inv. 


Bases 
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1 
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R 
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1 


1 


1 


1 


1 


1 


1 


-1 


-1 


a 


s z 


r 3 


1 


1 


-1 


1 


1 


-1 


-1 


1 


-1 


a 


zS z 


r 4 


1 


1 


-1 


1 


1 


-1 


-1 


-1 
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r 5 


2 


2 


-2 


-1 


-1 


1 


1 








a 


(S x -iS y ), 
-{S x -iSy) 


r 6 
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2 


-1 


-1 


-1 


-1 








a 


r 3 xr 5 


r 7 


2 


-2 





1 


-1 




V3 








c 


#3.1) 
























r 8 


2 


-2 





1 


-1 


V3 










c 


r 7 xr 3 


r 9 


2 


-2 





-2 


2 














c 


#1.-1). 
#1,1) 



We observe that: = {E; 2C 3 ; 3C2; cr h ; 2S* 3 ; 3d v } 

a. The order of Z) 3 h is 12. The representations l^V-Tg satisfy (5.49) by 
2 2 + 2 2 + 2 2 = 12. 

b. The elements a\ y and <7h belong to the same class because the C2 axes 
lie in the Uh plane. The "vertical" a v and a v belong to the same class 
both because of the C2 and because of a^. C2 and C2 belong to the 
same class because of either <7h or one a v . 

c. By Opechowski's theorem, C 3 and C3 belong to different classes. So 
do 5 3 and S" 3 . The element E forms its own class. This gives the three 
classes, imposed by k = r, eq. (5.54). 
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5.3 Construction of Relativistic MO:s 
5.3.1 Projection Operators 

The functions (f>\, belonging to row A of irrep k can be projected out from 
an arbitrary function by using the projection operator 

^ } = |E r S(^A (5.56) 



A 



Here g = jfcG, n^= dimension if irrep k, T^ k \A) is the representation matrix 
for group element A in the irrep k and Oa is the symmetry operation A. 

One can first operate on an arbitrary function, tp, by P^j obtaining 

4>1 = PgV (5-57) 
One then obtains the partners of <p\ by 

4>t = PPM (5-58) 

Suppose now that the basis functions in our LCAO-MO consist of atomic 
spinors, 

The bispinors Xjlm = \l\j m ) transform under the rotations according to 

As noted by Rosen and Ellis (1979), then the result of the projection 
operator becomes 

P^Xllm =|E ^(AY £ Xjlm'i-y^L'm^) (5-60) 
y A=\ m'=-j 

Here ta = 1 if A contains the inversion operator / and ta = otherwise. 
Note that, outside the symmetry center, the result of OAtp^i m may be shifted 
from the atom t to other, symmetry-equivalent atoms. 

Example 7: 

(Pyykko and Toivanen 1977): Relativistic symmetry orbitals for the 
double group D^. The origin is situated at the symmetry center. The 
variable s is +1 for k < and -1 for k > 0. 
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irrep even I odd / 
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Latest program: J. Meyer et al. (1996). 



5.3.2 Coupling Constant Method 

For the required coupling constants, see Koster et al. (1963) or Altmann 
and Herzig (1994). One now first constructs the non-relativistic MO:s, u\, 
and then couples them with the spin functions, 



0?i 



/2 



1/2 



1 

'2' 
1 

'2' 





a, 



(5.61) 



using the coupling constants or Clebsch- Gordan coefficients u l ^ v : 



ij,k 
Xfi,nu 



u\vl 



(5.62) 



The non-relativistic funcitons u\ are expressed in the \lmi) representa- 
tion. From that \lmi\m s ) representation, (5.62), one can pass to \l\jm,j) 
representation using the Clebsch-Gordan coefficients 



21+1 
l-m+\ 
21+1 



l-m+\ 
21+1 
l+m+\ 
21+1 



This method will automatically give the correct phase between the j = 
I ± I components. 
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5.4 Time Reversal 

5.4.1 Non-relativistic Case 

Time reversal is obtained by complex conjugating the wave function, 

Tip = ip* = Kip. (5.63) 
Let T = K (complex conj.), [H,T] = 0. 

(m-)V = #V ^ (m^y)^* = T(H1>) = H(Tip) = Hr (5.64) 

5.4.2 Inclusion of Spin 

Now 

T%P{r, s) = (-i) 2 >(r, -s)*. (5.65) 
This operator T is antilinear, 

r(ai^i + a 2 V> 2 ) = a^T^i + (5.66) 

and antiunitary, 

(TVi|TV>2> = (^IV'i)- (5-67) 

5.4.3 n-electron Wave Functions 

We now can write T as the product of (antiunitary) operator K and a 
unitary operator U. Possible choices are 

T = a ly a 2y .. .a ny K, a y = I ° 1 I (5.68) 



u 

or 



i 0, 



Tipiz^S!,. . . ,r n ,s n ) = (-i) (s ± + . . . + s n )ip(r l ,-s 1 ,. . . ,rn,-s n )*. (5.69) 
5.4.4 Kramers' Theorem 

For a system with an odd number of electrons, ip and Tip are orthogonal. 
In general, for arbitrary ip and <fi, 

(T^\T<f>) = (UK^j\UK(f)} = (KiP\KiP) = (5.70) 



<TV#> ( = 0) (TiP\T 2 iP) 



T 2 ip = -ip 



= (Tip\-ip) = -{Tip\ip) = o (5.71) 

Note that the equivalence used, T 2 ip = —ip, only holds for an odd number 
of electrons! 

In the absence of magnetic field, all energy levels of such systems are at 
least doubly degenerate. 
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5.4.5 The Cases (a), (b) and (c) of Wigner (1932) 

Consider the representation matrices, V, of a group G*. The operator T 
either does or does not give new, degenerate functions, spanning the repre- 
sentation matrices T*. The following possibilities exist: 

a. r may be real. 

b. £ is complex and non-equivalent to T* 

c. T is complex and equivalent to F* 



a. The degeneracy is doubled (from n for T to In. The functions ip and 
Tip belong to the same irrep. 

b. The degeneracy is doubled. If tp spans the irrep T, %p* spans another 
irrep Y* . 

c. The degeneracy remains n. 

According to the Frobenius-Schur criterion, 



Then 



2$, 



case (a) 
case (b) 
case (b) 



2n 








(5.72) 



A=l 



I -2g, 



5.4.6 Further Examples 



Some other examples of time inversion are considered below 



Momentum 




(5.73) 
(5.74) 
(5.75) 



Tp = —pT, the sign of p is inverted! 
TpT~ x = -p 



Angular momentum 



I = r x p | (5.74), p changes sign under time inversion 



Tl = -IT 



(5.76) 
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Spin- orbit Hamiltonian 



K 



hso = ^ 2 (— x ' P | both cj and p change sign 



5pm functions 



., 2 .i fl (5.69) ./0 



Ta = (— v = -i 



T/3 = (-ifH)« = ifi 



-i 

; i °/ 

-i 

1 0; 



1\ (5.68) 

0/ 



Rosch: T 



-1 

1 0, 



Eschrig: 



Ta = p 
Tf3 = -a 



1 
-1 

1 
-10/ 



-1 

-c(p x - \p y )/A 



i 
o 

cp z /A 

\c(p x + iPy)/A) 

\ 



g-ip-r _ _ e -ip-r 



-ip-r 



' o N 



cfe - ip y )/A 
\ -c Pz /A J 
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Compare the above with (3.45). Still, T 2 = -1 

External electromagnetic fields 
Invert the currents, 

3_^-b A^-A (5.82) 

5.5 Quaternions 

For complex numbers z = x + ly = (x, y) both addition: 

(x,y) + (X,Y) = {x + X,y + Y) (5.83) 

and multiplication: 

(x, y) (X, Y) = (xX - yY, xY + yX) (5.84) 
are defined. The multiplication is 

a. commutative, z\Zi = Z2Z\ 

b. associative, a(bc) = (ab)c 

c. norm- conserving, \z\zi\ = \zi\\z2\, \z\ 2 = x 2 + y 2 

Question: Can (a)-(c) be satisfied for more than two ordered compo- 
nents? 

Answer: Hamilton and Cayley (1843): For quaternions q = (a, b, c, d), 
with \q\ 2 = a 2 + b 2 + c 2 + d 2 , a non-commutative product with (c) can be 
found. 
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